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Abstract 

We investigate in detail a recently introduced “coherent averaging scheme” in terms of its useful¬ 
ness for achieving Heisenberg limited sensitivity in the measurement of different parameters. In the 
scheme, N quantum probes in a product state interact with a quantum bus. Instead of measuring 
the probes directly and then averaging as in classical averaging, one measures the quantum bus or 
the entire system and tries to estimate the parameters from these measurement results. Combining 
analytical results from perturbation theory and an exactly solvable dephasing model with numerical 
simulations, we draw a detailed picture of the scaling of the best achievable sensitivity with N, the 
dependence on the initial state, the interaction strength, the part of the system measured, and the 
parameter under investigation. 
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I. INTRODUCTION 


Averaging data is a common procednre for noise reduction in all quantitative sciences. 
One measures the noisy quantity N times, and then calculates the mean value of the N 
samples. Assuming that the useful signal part is the same for each run of the experiment, 
the random noise part averages out and leads to an improvement by a factor y/N of the 
signal-to-noise ratio (SNR). Instead of measuring the same sample N times, one may of 
course also measure N identically prepared samples in parallel, in which case we will think 
of them as “probes”. A lot of excitement has been generated by the realization that in 
principle one may improve upon the \/N factor by probes that are not independent, but in 
an entangled state. It was shown jTj that with such “ quantum enhanced measurements” 
the SNR can be improved by up to a factor N. Unfortunately, on the experimental side, 
decoherence issues have limited the quantum enhancement to very small values of N 
H]. For practical purposes it is therefore often more advantageous to stay with a classical 
protocol and increase N j5]. Since the decoherence problem is very difficult to solve, one 
should think about alternative ways of increasing the SNR through the use of quantum 
effects. One such idea is “coherent averaging”. The original scheme, hrst introduced in iiizi 
and named as such in j2], works in the following way: instead of measuring the N probes 
individually, one lets them interact coherently with a iV + 1st system (a “ quantum bus”) and 
then reads out the latter. In this way, quantum mechanical phase information from the N 
probes can accumulate in the quantum bus, and this can improve the SNR also by a factor 
N, even when using an initial product state (see FigQ. A physical example considered in 
detail was the coupling of N atoms to a single leaky cavity mode, which allowed to measure 
the length of the cavity with a precision scaling as 1/N, which corresponds to the above 
SNR oc N. This scaling is the long-sought Heisenberg-limit (HL), contrasting with the 
I/a/IV scaling characteristic of the standard, classical averaging regime, also called standard 
quantum limit (SQL). 


So far, however, the method was limited to estimating a parameter linked to the inter¬ 
action of the N probes with the quantum bus. This makes comparison of the performance 
with and without the coupling to the quantum bus impossible, as in the latter case the 
parameter to be estimated does not even exist. In the present work we go several steps 
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Figure 1. Classical averaging (left) versus coherent averaging (right). In coherent averaging, the N 
probes are not read out individually and the results averaged, but one lets the probes coherently 
interact with a quantum bus, and then either measures the latter, or a global observable of the 
entire system. The parameter to be estimated can parametrize the probes, the quantum bus, or 
the interaction. 


further. Firstly, we extend the scheme to estimating a parameter that characterizes the 
probes themselves, or the quantum bus itself. Secondly, we analyze in detail conditions for 
the observation of the HL scaling by systematically studying strong, intermediate and weak 
coupling regimes. Numerical simulations are used in order to verify and extend results from 
analytical perturbation-theoretical calculations. Thirdly, we investigate the question which 
part of the system should be measured. 


Note that achieving HL scaling of the sensitivity in coherent averaging with an initial 
product state is not in contradiction with the well-known no-go-theorem [1] which is at the 
base of the often held believe that entanglement is necessary for surpassing the SQL. The 
reason is that in [Tj the Hamiltonian is assumed to be simply a sum of Hamiltonians of 
independent subsystems with no interactions, which is a natural assumption when coming 
from classical averaging. Meanwhile, however, several other ways have been found to bypass 
the requirements of the theorem and thus avoid the use of entanglement for HL sensitivity, 
notably the use of interactions (also known as non-linear scheme) |9l |T0], multi-pass schemes 
|1], or the coding of a parameter other than through unitary evolution (e.g. thermodynamic 
parameters such as the chemical potential) |ll| . 

From a perspective of complex quantum systems, the models that we study are typical 
decoherence models: the quantum bus may be considered an environment for the N probes, 
or vice versa. However, in general we will assume that we can control both probes and 
quantum bus, and in particular prepare them in well defined initial states which we take as 
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pure product states or thermal states. 


II. MODELS AND METHODOLOGY 


A. Models 


The systems we are interested in have the following general structure depicted in Figj^ 
The corresponding Hamiltonian can be written as 


H 


6Hq + eHi^t 


( 1 ) 

( 2 ) 


\ ^=1 / \ / 
where i^o contains the “free” part (probes and quantum bus), and ifmt the interaction be¬ 
tween the probes and the quantum bus. We have introduced two dimensionless parameters 
6 and £ which we will use to reach the different regimes of strong, intermediate, and weak 
interaction. In the second line we specify the Hamiltonians for N non-interacting probes 
which we assume to depend on the parameter cui, and the Hamiltonian Hr of the quantum 
bus (or “reservoir” in the language of decoherence theory) which depends on the parameter 
cjo- The interaction has the most general form of a sum of tensor products of probe-operators 
and quantum-bus-operators and we assume that it depends on a single parameter x. 


As specific examples of systems of this type we consider spin-systems, where both the 
probes and the quantum bus are spins-1/2 (or qubits) and thus described by Pauli-matrices 
A, Y, Z. Without restriction of generality, we can take Hi = for the i-th probe, and 

Hr = (h = 1 throughout the paper), where the bracketed superscripts denote the 

subsystem, and the zeroth subsystem is the quantum bus. For the interaction we focus on 
two different cases: an exactly solvable pure dephasing model with 

= (3) 

i 

and a model that allows exchange of energy through an XX-interaction, given by 

Hi,t = |5^X«0XW. (4) 

i 

We refer to these two models as ZZZZ and ZZXX models. 
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B. Initial state 


Given the difficulty of producing entangled states and maintaining them entangled, we 
consider here pure initial product states with all the probes in the same state, which may 
be different from the state of the quantum bus. For the spin-systems we parametrize these 
states as 


I V’o) 


I ® 10 

(cos(q!)| 0 ) -I- sin(Q!)e‘'^| 1 
0 (cos(/9)| 0 )sin(/9)e“^| 1)) , 


( 5 ) 

( 6 ) 


where | 0), 11) denote “computational basis states”, i.e. 2’| 0 ) = | 0) and 'Z\\) = — | 1) for 
any spin. Eq.(|^ implies that in the subspace of the probes, the initial state is a SU{2) an¬ 
gular momentum coherent state of spin j = N/2. Since both initial state and the considered 
Hamiltonians are symmetric under exchange of the N probes, this symmetry is conserved 
at all times, and allows for a tremendous reduction of the dimension of the relevant Hilbert 
space: from 2^^^ to only 2{N -|- 1) = 2{2j + 1) dimensions. The corresponding basis in the 
probe-Hilbert space is the usual joint-eigenbasis \j,m) of total spin and its 2 ;-component. 
We will omit the label j and have thus the representation of | "00) in the symmetric sector 
of Hilbert space 

I V’o) = 

(cos(/?)| m, 0)-I-sin(/3)e'‘^| m, 1)) . (7) 



For the ZZZZ model we also consider thermal states of the probes, see eq.(37) below. In other 
contexts, the above models have been called spin-star models, and analyzed with respect to 
degradation of channel capacities and entanglement dynamics [UMI. 


C. Quantum parameter estimation theory 

The question of how precisely one can measure the parameters uJi,uo and x is addressed 
most suitably in the framework of quantum parameter estimation theory (q-pet). Q-pet 
builds on classical parameter estimation theory, which was developed in statistical analysis 
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almost a century ago [T^ [T6] . There one considers a parameter-dependent probability dis¬ 
tribution p{A, 6) of some random variable A. The form of p{A, 9) is known, and the task 
is to provide the best possible estimate of the parameter 9 from a sample of n values Ai 
drawn from the distribution. For this purpose, one compares different estimators, i.e. func¬ 
tions 6'est(^i, • • • ,An) that depend on the measured values Ai (and nothing else), and give 
as output an estimate 9est of the true value of 9. Since the Ai are random, so is the estimate. 
Under “best estimate” one commonly understands an estimate that fluctuates as little as 
possible, while being unbiased at the same time. 

In quantum mechanics (QM), the task is to estimate a parameter 9 that is coded quite gen¬ 
erally in a density matrix, p{9). One has then the additional degree of freedom to measure 
whatever observable (or more generally: positive-operator valued measure (POVM)|17)). 
The so-called quantum Cramer-Rao bound is optimized over all possible POVM measure¬ 
ments and data analysis schemes in the sense of unbiased estimators. It gives the smallest 
possible uncertainty of ^est no matter what one measures (as long as one uses a POVM 
measurement — in particular, post selection is not covered, see [12] for an example), and no 
matter how one analyzes the data (as long as one uses an unbiased estimator). At the same 
time it can be reached at least in principle in the limit of a large number of measurements. 
The quantum Cramer-Rao bound (QCR) has therefore become the standard tool in the held 
of precision measurement. It is given by 

Var(0est) > ^ , (8) 

where Var(6'est) is the variance of the estimator, Ig the Quantum Fisher Information (QFI), 
and M the number of independent measurements. A basis-independent form of Ig reads 


Ig — 2 


dstT[dgpge P'^^'dgpge 


In the eigenbasis of pg, i.e. for pg = YlirPr \ V’r )(V’r | we obtain 

{dgPrf , ^'^{Pn- Pmf 








( 9 ) 


( 10 ) 


where the sums are over all r and n, m such that the denominators do not vanish. It 
is possible to give a geometrical interpretation to the QFI, namely in terms of statistical 
distance. To this end one dehnes the Bures distance between two states p and a as 


ds(p, a) = \/2^ 1 - tr[(pV2^pi/2^i/2j^ _ 


( 11 ) 
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In the case of two pure states 0, V’, we have 


0), I V')) = ^2^1- \ {(j)\'ip)\) . (12) 

The Bures distance was shown to be related to the QFI by |2n| 

ie = Adlipie),pie + de))/ 4 . ( 13 ) 

It provides an intuitive interpretation to the best sensitivity with which a parameter can 
be estimated in the sense that what matters is how much two states distinguished by an 
infinitesimal difference in the parameter 6 differ, where the difference is measured by their 
Bures distance. In the case of a pure state, the QFI is equal to 

/, = 4(( de m \de'iP{9))-\{ 'iP{9) \de'iP{9))\^) . (14) 

D. Perturbation theory 

It is clear that the model ([^ cannot be solved in all generality. One way of making 
progress is to use perturbation theory. This can be done in two ways: In the standard use 
of perturbation theory one solves the Schrodinger equation for the free Hamiltonian Hq and 
then treats the interaction Ffint in perturbation theory, provided that the interaction is small 
enough. In the regime of strong interaction, one can do the opposite thing: solve the pure 
interaction problem first, and then calculate the additional effect of the free Hamiltonian 
as a perturbation. Formally this does not make a big difference. More important, already 
on the level of the expression for the QFI, is the question whether the parameter to be 
estimated enters in the perturbation or in the dominant part of the Hamiltonian. We call 
the perturbation theory relevant for these two cases PTl and PT2, respectively. 

To better understand the difference between PTl and PT2, consider a Hamiltonian con¬ 
taining two parts where one of them depends on a parameter that we want to estimate, 

H{9) = H^{9) + H2, (15) 

and the state 

l^(^)) = exp(-hiP(0))| V'o) • (16) 
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In PTl we switch to the interaction picture with respect to ip 2 , 




( 17 ) 


Under the conditions that ^ 1 and ||i7i(6')|| -C ||-f72|| we can use second 

order perturbation theory in order to calculate the QFI 

ie = ^[ [ dhdt2K\^,){H,y{e,h),H,y{e,h)), (is) 

Jo Jo 

with K\^){A, B) = {yi AB \ y) — {y\A\y){y\B\y). if the free Hamiltonian commutes with 
the interaction Hamiltonian (as is the case in the ZZZZ model), and under the assumptions 
that [Hi{6), Hi{6)'] = 0, one can calculate the QFI exactly, 

ie = 4 yKi^,){Hy{e),Hy{e)). ( 19 ) 


The last term of the right hand side is the variance of Hi (6) in the initial state, and we 


thus recover a well known result in Q-pet |2T]. At the same time, eq.(18) tells us that the 
QFI is of second order in the perturbation. 


In PT2, we do the opposite: estimate the parameter linked to the Hamiltonian that 


dominates. With the notation of eq.(15) this means that we try to estimate the parameter 


0, considering that \H 2 jt\, \H 2 yt\ 1 and ||i7i(6*)|| S> ||772||, with 




( 20 ) 


In this case, the result (18) does not apply anymore. Indeed, to obtain (18) we calculated 


the overlap {y{9, t) \ y{6 + dO, t )), which equals ( yij{9, t) \ yij{9 + d9, t)). Now we need 


{ t) |e '”0«)e + ii0, t )). 


( 21 ) 


Here we have dehned 


V'i,/(^U)) 


T 



Hj,i{t')dt') 


I V’o) 


L -'0 J 

for j = 1, 2, and T is the time-ordering operator. In this case, the lowest order term that 
appears in the expansion of the QFI is the unperturbed term. 


h = ie {{,p„ I HP(ef I ,p„} - I HP(B) I v>o>"). 


( 22 ) 








The formal range of validity is now given by \H 2 jt\,\H 2 j't\ 1 and ||i/i(6*)|| S> ||-ff 2 ||- 

The first and second order terms are too cnmbersome to be reported here. Bnt since in the 
formal range of validity of the perturbation theory they have to remain small in comparison 
to the zeroth order term, the scaling of the QFI with N will be given anyhow by the zeroth 


order term (22). 

For practical applications and in the spirit of the original “coherent averaging” scheme, 
we are also interested in the sensitivity that can be achieved by only measuring the quantum 
bus. To this end, we calculate the reduced density matrix by tracing out the probes, and then 
the QFI for the corresponding mixed state which we call “local QFI” lf'\ Since tracing out a 
subsystem corresponds to a quantum channel under which the Bures distance is contractive 
[22], we have h > lf'\ 

In general, the calculation of the QFI for a mixed state is rather difficult, as one has to 
diagonalize the density matrix twice (either for two slightly different values of the parameter 


for calculating the derivatives (10), or for calculating and Techniques 

j. Here, 


for bounding the QFI for mixed states have been developed in the literature 
we only calculate the QFI for the mixed state of a single qubit, which is easily achievable 
numerically. 

Another important practical question is with which measurement the optimal sensitivity 
can be achieved. In principle, the answer can be found from the QFI formalism [I9[ 1^ . 
Here we follow the strategy of considering local von Neumann measurements of the quan¬ 
tum bus and comparing the achievable sensitivity to the optimal one. If a von Neumann 
measurement with a Hermitian operator B is performed, one can show that an estimator 
6*est(5i,..., Hm) = Bi/M) with f{6) = {i/(6)\El'll;(6)) leads to first order in the 

expansion of f~^ to an uncertainty (standard deviation) of ^est given by 


' v'M\de{B}\ 


with Var(i?) = ( H (see eq.(9.16) in |25]). This “method of the first moment” can 

always be rendered locally unbiased by adding a shift to 6*est- The uncertainty corresponds 
to the minimal change of 9 that shifts the distribution of the by one standard deviation, 
assuming that the shift is linear in dO for small dO. Since df is based on a particular 
estimator, we have < MIq. 

We call the local observable of the quantum bus = 11®'^ ® A, and and the 
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corresponding uncertainties for the estimation of x and oji. The QFI of the reduced density 
matrix for the quantum bus alone will be denoted and We have MIq > > 

(5^)“^ = The last step follows from = (A) for any quantum state. 


E. Numerics 

When ehfint and 5Hq are of the same order, both forms of perturbation theory typically 
break down. Unless one has an exactly solvable model (such as the ZZZZ model), one has 
to rely on numerics. In addition, we use numerics to test all our analytical results. The 
perturbative results are, in general, limited to a hnite range of iV, such that when one 
wants to make a statement about the scaling of the sensitivity of a measurement with N for 
large iV, one has to rely once more on numerics. All numerical calculations use one of the 
spin-Hamiltonians, eq.([^ or (|^. 


The numerical results are obtained by calculating the time evolution operator U{t) = 
exp(—ihft) for the full Hamiltonian in the Schrodinger picture, propagating the initial state 
([^ for two slightly different values of the parameter we are interested in (a:, wi, or wq), 
obtain from this a numerical approximation of the derivatives of |'0(^)) with respect to the 


parameter, and then calculate the overlaps in eq.(14). In this way, we obtain the “global 
QFI”, which is relevant if one has access to the entire system (i.e. probes and quantum bus). 
To check the stability of the numerical derivative, we calculate numerical approximations of 
the derivative for two different changes in the value of the parameter, 10“® and 10“®. 

For the spin-Hamiltonians considered, the reduced density matrix of the quantum bus 
is the density matrix of a single spin-1/2 which simplihes the calculation of the QFI. For 
numerical calculations we use the basis-independent form ([^ of the QFI and perform the 
integral analytically. We also calculated and 5^^°^ numerically “exactly” by directly 


evaluating (23). 


In order to check the validity of the perturbative result, we verihed that in the range of 
validity of the perturbation theory the difference between the exact QFI and the perturbative 
result scales as 5^ or as function of the perturbative parameter 6 or e. 
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III. RESULTS 


We now present our results for the estimation of x, ooo, and oji in the different regimes, 
focusing hrst on the global QFI. All hgures shown have the parameters uq = 1, wi = 1, 
f = 1, X = 1. The initial pure state ([^ is taken always with a = tt/S, /? = vr/h, 0 = Stt/S, 
= Stt/S unless otherwise indicated. 


A. Global QFI 


1. Estimation of x 

The perturbation theory for estimating x for small interactions was developed in 
Inserting the form (|^ in eq.(18), one hnds that for identical and identically prepared systems 


<Si = Si, and \ip)i = | )) the correlation function in eq.(18) is given to lowest order in 

e by 

■^1 ho ) ( 3 ^) U)) (®)^2)) 

{SJix, ti), S'/(x, fa)) (C I Ru{ti)Rf,{t2) \ 0 




+ N‘^{ip\Si,'{x,ti)\ip){ip\Sf,'{x,t2)\y^)K\oiRAti),Rni't‘i))} ■ (24) 

We have dehned TTmt,/ = Uo{t)eHi.aRJo{ty with U^it) = exp(i5hfo^)) and S^{xR) = 
Uo{t)S^UQ{t)\ R^{t) = UQ{t)R^UQ{ty. This implies a structure R = e^{ni^xN + n 2 ,xN‘^) + 
0{e^) of the QFI, where ni^x and n 2 ^x can be expressed in terms of time integrals of corre¬ 
lation functions. However, higher orders in e: limit the formal validity of the perturbation 
theory to sufficiently small values of N. Indeed, the next higher order may contain terms of 
the order e^N^, which are only much smaller than the second order for N <C Xje. 



0 JO 


Eq.(24) allows one to establish the condition for HL scaling, namely that (7] 

dt\dt2 ^((p I Sj{x, h) I (p)((p I V(u ^ 2 ) I 

fl,U 

xAig) {Ru{ti), Rij,{t2)) 7 ^ 0 . 

Numerics for the ZZXX model con&ms the perturbative result in its expected range of 
validity. Moreover, it also indicates that the HL scaling works beyond the formal range of 
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Figure 2. From left to right: Global QFI for the ZZXX model for x for weak, medium, and strong 
interactions (e = 0.001, 1, e = 100, and 6 = 1). Blue X-symbols: exact numerical results. Purple 
circles: perturbative result (PTl). Red diamonds: zeroth order (unperturbed) term in PT2. The 
dashed orange (resp. green continuous) lines represent f{N) oc N'^ (resp. N). 

validity of the PT. This is shown in Figj^ where we compare the global QFI for measuring 
X for weak, medium, and strong interactions. We see that PTl works correctly for eN -C 1. 
For medium and strong interactions, PTl still predicts a scaling of the global QFI propor¬ 
tional to N"^. While this is conhrmed by the exact numerical results, the prefactors differ 
outside the formal range of validity of perturbation theory. The scaling is more easily 
observed for strong interactions than for weak ones, but Figshows that even for weak and 
medium interactions a N'^ component is already present. The onset of this behavior can 
clearly be identihed in Figj^for e = 0.001 and e = 1. 


For strong interactions, PT2 is appropriate for obtaining the QFI for x. The zeroth term 


(22) dominates in the range of validity of the perturbation theory, and leads to 


Ix = (iV^ sin^(2Q;) cos^(2/3) -|- N cos^(2q;)) -|- 0{5). 


(26) 


This implies HL scaling in the formal range of validity {N5 <C 1 and 5 e). Figure]^ 


shows that eq.(25) works well even beyond this regime. A more precise assessment of the 


range of validity has to consider the matrix norm of H^t. If the largest absolute eigenvalue 
of Hi is Amax, then PT2 is expected to work for <C 1 and 5 e. 


2. Estimation of oji 

The situation is similar for estimating ui. PTl (i.e. treating 6Ho as perturbation, such 
that in the interaction picture Hq j = f/int(t)(5hfofAint(Q^ with f/int(Q = e‘^*^‘"‘), and assuming 
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Figure 3. From left to right: Global QFI for the ZZXX model for oji for weak, medium, and strong 
interactions (h = 100,1, <5 = 0.001, and e = 1) . Blue X-symbols: exact numerical results. Purple 
circles: perturbative result (PTl). Red diamonds: zeroth (unperturbed) term in PT2. The dashed 
orange (resp. green continuous) lines represent f{N) oc (resp. N). 


that [Ru, = 0, Vz/, /i, leads to a correlation function to lowest order in 5 given by 




( 26 ) 




( 0 ), 


with [ujiRi) = f/int(ti)hf/(c(;i)t/j^^(ti). Note that is still an operator on the 

quantum bus after sandwiching it between probe states \^)- Eq.(26) together with (18) 
shows that obeys HL scaling for 5 N l.It is easier to observe HL scaling for 5 -C 1, 
i.e. in the regime of small free Hamiltonian or, equivalently, strong interactions, see Figj^ 
For medium interactions (5 = 1), HL scaling is still observed, whereas for weak interactions 
(5 = 100) SQL scaling prevails at least up to iV = 2000. Formally, the range of validity of 
PTl is limited here to iV 1/5, but numerics indicates HL scaling up to much larger N. 

In the formal range of validity, PTl gives the necessary condition Hi] 7 ^ 0 for 
observing HL scaling, as otherwise becomes proportional to the identity operator in 
the quantum bus Hilbert space. Numerically it can be checked that a violation of this 
condition indeed leads only to SQL scaling. We verihed this for the ZZZX model, dehned 
as the ZZXX model, but with a Hamiltonian Si = |Zk) instead of Si = and with 

5 = 0.001, 0.1,1 and 5 = 100 for several random initial states. 

For weak interactions, using PT2, one finds a QFI with the structure = aoe^N + 
aie^N'^ + a 2 S^N^ + 0{e^) with some coefficients a^. In the range of formal validity, the 
dominating term is e^N. This once more implies that SQL scaling dominates the estimation 
of Ui for weak interactions, in agreement with the left plot in Figj^ 
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Figure 4. From left to right: Global QFI for the ZZXX model for ojq for weak, medium, and strong 
interactions {5 = 100, 1, h = 0.001, and e = 1). Blue X-symbols: exact global QFI for wq. The 
dashed orange (resp. green continuous) lines represent f{N) oc N'^ (resp. N). Same state as in 
FigjU 


3. Estimation of uiq 

In Figj^we show numerical results for for the ZZXX model. We see that for uq the 
coupling of additional qubits to the central one not only does not improve the best possible 
sensitivity for N larger than a number of order one, but in general even deteriorates it. A 
perturbative analysis in the framework of PTl is not helpful here, as (interaction 

picture with respect to is an operator that acts non-trivially in the full Hilbert space. 


B. Local QFI and local quantum bus observable 

As we have seen in the last section, the global QFI indicates that HL scaling can be 
observed with a Hamiltonian of the form ([^ and an initial product state. We now investigate 
whether it is enough for achieving the HL to measure the quantum bus. To this end, 
we calculate the QFI of the reduced density matrix of the quantum bus, as well as the 
uncertainties of the parameter estimates based on a specihc observable A of the quantum 
bus. We do not investigate further the estimation of uo as already the global QFI shows 
that the sensitivity cannot be improved by coupling to additional qubits. 
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N N N 

Figure 5. Local QFI and inverse squared uncertainties of x based on the local observable = 
(^(0) _|_ Z^^'^)/2) for the ZZXX model for weak, medium, and strong interactions, £ = 0.001, 0.1, 
£ = 100, and <5 = 1 from left to right; Blue X-symbols: exact numerical result for Purple 

circles: perturbative solution for Red crosses: exact solution for The dashed 

orange (resp. green continuous) lines represent f{N) oc N‘^ (resp. N). Same state as in Figj^ 


1. Estimation of x 


The behavior of 5^^°^ was analyzed in second order perturbation theory in [7j. Within 
its range of validity, HL scaling was found under the condition of a noiseless observable of 
the quantum bus that remains noiseless without interaction with the probes. Here we relax 
the conditions and give a more general form in the appendix, eqs.( 40pl| ) together with (23). 
Figj^ shows HL scaling of the sensitivity for weak interactions (e: = 0.001) for N < 500 and 
the measurement of the quantum bus + Z^^'>)/2). The perturbative result for 


( 5 : 


A ( o )\-2 


agrees perfectly well with the exact numerical result in this regime. We also see 


that the local QFI provides an upper bound to However, the local QFI is rather 

small in the range of N accessible to exact numerical evaluation, such that at least for these 
values of N the observed HL scaling is not of much use. For larger e, PTl quickly breaks 
down, as is shown in Figj^for medium and strong interactions: For £ = 0.1 the break down 
occurs at ~ 10, compatible with Ne ~ 10. For e = 100, PTl is already invalid in the 
sense that the QFI becomes negative at = 1 and we do not plot it. Moreover, the exact 
numerical values both for and show that for strong and medium interactions the 

achievable sensitivity through the measurement of the quantum bus alone deteriorates with 
increasing N for large enough N. 
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Figure 6. Local QFI and inverse squared uncertainties of coi based on the local observable = 
for the ZZXX model for weak, medium, and strong interactions, 6 = 100, 1, 
6 = 0.001, and e = 1 from left to right. Blue X-symbols: exact numerical result for Purple 
circles: perturbative result for Red crosses: exact result for The dashed orange 

(resp. green continuous) lines represent f{N) oc N'^ (resp. N). The saturation at 10“^® reached in 
the right plot corresponds to the numerical precision. Same state as in Figj^ 


2. Estimation of oji 

A perturbative result for could be obtained in the case = 0, [Hn^Ry] = 

0,Vz/, and [Ry,R^ = 0,Vi/,/x. The hrst condition avoids that in the interaction picture 
acts non-trivially in the full Hilbert space. The second and third conditions avoid that Hq 
in the interaction picture acts non-trivially in the full Hilbert space. However, all three 
assumptions together lead to a diverging as they imply = 0. Numerically we 

can explore a more general case where these conditions are relaxed. The results are shown 
in Fig. For strong and medium interaction (5 = 0.001 and 5 = 1), while the global QFI 
shows HL scaling, the local QFI goes to zero, leading to the impossibility to estimate cui by 
a local measurement. For weak interactions (5 = 100), the local QFI shows neither clear HL 
nor SQL scaling. Nevertheless, the SQL behavior of the global QFI sets an upper bound on 
the local QFI. 


C. Exact results for the ZZZZ model 

1. Pure product state 

In order to corroborate the above results we calculated the QFI for the different parame¬ 
ters and and exactly for the ZZZZ model. The expressions of the global QFI 
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for the initial state ([^ are given by 

Ix = cos^(2q!) sin^(2/?) + Nt'^e^ sin^{2a) 

luji = sin^ {2a) 

Lq = 5H'^siv?{2l3). 


(27) 

(28) 
(29) 


Ix clearly shows HL scaling as long as cos(2q!) sin(2/?) ^ 0, while cui can only be measnred 
with a sensitivity that scales as the SQL. The best possible estimation of ujq does not profit 
from the quantum probes at all as is independent of N and as for the ZZXX model we 
do not investigate it any further. All global QFIs show a scaling oc demonstrating that 
the sensitivity per square root of Hertz can still be improved by measuring for longer times, 
in contrast to the typical time dependence of classical averaging. 

The general results for the local quantities are cumbersome with the exception of 
which vanishes for all initial states (|^ as the reduced density matrix of the quantum bus 


does not depend on ui, see eq.(44) in the Appendix. For the estimation of x we give the 
reduced density matrix and the uncertainty obtained via a measurement of in the 
appendix. Here, we provide results for two specific initial states. The most favorable case 
for the estimation of x, a = 0, (3 = 7r/4, i.e. 


|V’o) = |O)®'^0(|O) + |l))/\/2, (30) 

leads to the global QFI Ix = For the local QFI we have I^'* = We 

notice that Ix = l‘x \ i-e. restricting ourselves to a measurement of the quantum bus does 
not affect the best possible sensitivity for a estimation of x, and that precision follows HL 
scaling. Moreover, one can easily show that the corresponding QCR bound is reachable by 
measuring X(o). 

Now consider the initial state with a = 7r/4 and 13 = 7r/4, and (j) = (p = 0, i.e. 

I ^0 ) = (I 0 ) + 11 ))®^ ® (I 0 ) + 11)) /2(^+')/2 _ (3^) 


This is the worst pure state for measuring x. We obtain 


I„ = Nt^E^ , 


(32) 


for the global QFI, i.e. SQL scaling, and 


NH'^e^ isji^{etx) 
cos{etx)~^^ — 1 ’ 


(33) 
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for the QFI of the quantum bus. For etx = 71/2 the QFI vanishes, which can be understood 
from the fact that the reduced density matrix does not depend on x. For the uncertainty of 
X based on the measurement of we hnd the exact result 


exact 


-2 


tax? {etx) 


cos{etx)~'^^ cos{6ooot)~‘^ 


(34) 


This shows that for the initial state (31), both the local QFI (33) and ((5,^exact)~^ de¬ 
cay exponentially with N for sufficiently large N, i.e. this state is not suited for coherent 
averaging if we can only measure the quantum bus. 

s:X(0) 


It is instructive to use PTl for calculating (5,Pexact) which leads to 

NH^e^x^ 


^-2 _ 
v^x,pert/ 


N?x'^e‘^ + tan^(5ci;oQ 


(35) 


If we expand in powers of e, we hnd + C^(^"^)- The exact result for 

(^^exact)”^ rehects the behavior of ji;°\ whereas the perturbative version, ((5,^pert)~^) predicts 
a completely different result, namely a scaling oc N for large N. If one stays in the range of 
validity of PTl, one does not notice that the uncertainty diverges. 


Therefore, for the initial state (31), the validity of the perturbative expressions for the 
uncertainty of an observable of the quantum bus and the local QFI does break down outside 
the range of validity of PTl, in contrast to the global QFI, where the perturbative expression 
still predicts the correct scaling behavior and only differs in the prefactor from the exact 
result. The decaying local QFI shows that the coherent averaging scheme does not allow 
one to reach HL scaling for the estimation of x through a measurement of the quantum bus 
only. 


In order to hnd out how generic the decaying local QFI is for different initial states, we 
investigated the dependence of the scaling on N on the angle a that dehnes the state of 
the probes, eq.([^. We keep ? = 7r/4, 0 = (^ = 0. Figure shows the local QFI for a: as a 
function of a and N. We see that when increasing a from zero, the QFI starts to decrease 
with N for N larger than some bound iVo(Q!), like in the case just studied {a = vr/4). The 
hgure also indicates that with increasing N the range of a leading to a non-decreasing QFI 
is reduced more and more. This shows that over an ensemble of initial states, a local QFI for 
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Figure 7. Local QFI for the ZZZZ model for x as a function of a and N, with e = S = 1, 
t = 1, X = 1, Wo = 1, = 1, /3 = vr/4, (/) = 0, (/? = 0. The contours are at = 

{ 10 ^^, 10 ^, 1 , 10 -^ 10 "^°, 10 - 30 , 10 - 60 , 10 - 166 , 10 - 160 }. 

the estimation of x that decreases with N is the norm, and the HL scaling for the optimal 
state an exception. 


S. Thermal state for the probes 


In order to answer the question how a lack of purity of the initial state affects our results, 
we take the N probes in a thermal state 


Pth — 


1 

Z 




(36) 


with Z = = l/^ksT), where T is the temperature and ks the Boltz¬ 

mann constant. The quantum bus is in a pure state | 'i/'bus) = cos(/9)| 0) -|-sin(/3)e 1), and 
the new initial state is the mixed product state 


do Pth ® I V^bus ) ( t/^bus 


(37) 


This resembles the DQCl protocol in quantum information processing that starts with all 
qubits but one in a fully mixed state, but which still allows one to solve a certain task more 
efficiently than with a classical computer (the “power of one qubit”) ^\. 
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The exact results for the global QFI read 

Ix = sin^(2/3)£^t^(A^^ tanh^(/9thi^i) + iV(l — tanh^(/3thi^i))) 

= A^/3th(l - tanh^(/3thWi)), and (38) 

I LOO = ^^t^sin^(2/3). 

This shows that it is possible to reach the HL scaling for the estimation of x using thermal 
states of the probes, even though the prefactor of the term becomes small for large 
temperatures (/Sthi^i •Cl). The level spacing of the probes can only be estimated with a 
sensitivity scaling as the SQL, and the thermal probes are entirely useless for improving the 
estimation of the level spacing of the quantum bus. 

Remarkably, the reduced density matrix has the same form as the one for the pure product 
state ([^ when setting 

cos^(q;) = e and sin^(a)^ = (39) 


This implies that for any pure product state ([^ there exists a thermal state of the probes 
with the same and hence the same best possible sensitivity of estimating x by measuring 
the quantum bus. For locally estimating cui, a thermal state of the probes is advantageous 


compared to the pure state (30) where the corresponding local QFI vanishes. The thermal 
state introduces a dependence on ui through the initial state that is absent for the pure 
states considered. If also the quantum bus is in a thermal state initially, the interaction 
strength cannot be measured in the ZZZZ model. 


IV. SUMMARY 

In summary, we have examined in detail a coherent averaging scheme for its usefulness of 
Heisenberg-limited precision measurements. In the scheme, N probes that are initially in a 
product state, interact with a quantum bus and one measures the latter or the entire system. 
Combining analytical results from perturbation theory and an exactly solvable dephasing 
model with numerical results, we have shown that this setup allows one to measure the 
interaction strength and the level spacing of the probes with HL sensitivity if one has access 
to the entire system. Strong interactions favor better sensitivities in this case. If one has 
only access to the quantum bus, the results depend on the initial state, but HL sensitivity 
is achievable only for the interaction strength and a small set of initial states. 
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Remarkably, for measuring the interaction strength in the exactly solvable ZZZZ model, 
there is a mapping of the local quantum Fisher information for thermal states of the probes 
to the one for pure states. Globally HL sensitivity for estimating the interaction strength 
can be achieved with thermal probes at any hnite temperature, as long as the quantum bus 
can be brought into an initially pure state. The sensitivity of measurements of the level 
spacing of the quantum bus cannot be improved by coupling it to many probes, even with 
access to the entire system, and in fact deteriorates with an increasing number of probes. 
Altogether, our investigations have led to a broader and more detailed view of the usefulness 
of the coherent averaging scheme and may open the path to experimental implementation. 


V. APPENDICES 

A. Uncertainty of a local observable 


If one relaxes the condition used in [ 7 ], namely that A| ^) = a^\^) and Hji] = 0, 
one hnds for the variance of the observable 


Var(A(°)) = (A2)-(A)2 + R / dhN{S,{h)){[RM, B]) 

Jo 


+e" 



dtidt2 {{N{N - l){S,y{ti) ){Si,{t2 )) 


0 JO 


+N{ S,{h)S,{t2) ))( [RM, B]R^{t2)) 

+{N{N - 1 )() ( S,{t2)) + N{ S,{t2)S,{h ))) 

{R^{t2)[B,R,{h)])} + e^ J J dhdt2N^{SM) 

{S,{t2)){[RAti),A]){[Rj2),A]) 


(40) 


where B = A^ — 2(A)A, the expectation values for A, B, and Rfj.{t) are taken with respect 
to I .^), and the expectation value for S^{t) = S^{x, t) is with respect to | 93 ). The derivative 
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of the mean value of is given by 


^ (ie j‘ dhN{S,(h)){lR^(,t,),A]} 

+ E^ f r dt,dh{{N(N - l){S,{t,)){S^{t2)) 
Jo Jo 

+ Su{ti)S^{t2) ))( [Riy{ti), A]Rfj^{t2)) 


+ {N{N - l){S,{h)) ( S^{t2)) + N{ S,{t2)S,{h ))) 

(fl„(fc)[4fl.(*i)]>}). 


From these two quantities we obtain 6 


A( 0 ) 

X 


according to eq. (23). 


( 41 ) 


B. Local analysis of ZZZZ 

The reduced density matrix for the ZZZZ model starting in a pure product state ([^ 
has the matrix elements 

Pod = cos^iP) (42) 

pdd = sin^(/?) (43) 

pl,°d = ^ sin(2/3)e + sin(a)2e ‘"^*)^ , (44) 

from which one can easily compute the local QFI. 

The relative uncertainty for x using a measurement is: 


('rX(O) \2 
\^x,pert/ 


1 - (^sin(2/3) Y11=-n/2 {m-N/ 2 ) cos(Q;)^+2msin('Q,)iv + P + 2extm)^ 

2etsm{2/3) Ylm=-N /2 cos(q:)^+2 ™ sin(Q;)^-2"^ sin((5a;ot + (p + 2extm) 

(45) 
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